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ABSTRACT 

Context. Earlier work has suggested that large-scale dynamos can reach and maintain equipartition field strengths on a dynamical 
time scale only if magnetic helicity of the fluctuating field can be shed from the domain through open boundaries. 
Aims. Our aim is to test this scenario in convection-driven dynamos by comparing results for open and closed boundary conditions. 
Methods. Three-dimensional numerical simulations of turbulent compressible convection with shear and rotation are used to study 
the effects of boundary conditions on the excitation and saturation of large-scale dynamos. Open (vertical-field) and closed (perfect- 
conductor) boundary conditions are used for the magnetic field. The shear flow is such that the contours of shear are vertical, crossing 
the outer surface, and are thus ideally suited for driving a shear-induced magnetic helicity flux. 

Results. We find that for given shear and rotation rate, the growth rate of the magnetic field is larger if open boundary conditions are 
used. The growth rate first increases for small magnetic Reynolds number, Rm, but then levels off at an approximately constant value 
for intermediate values of Rm. For large enough Rm, a small-scale dynamo is excited and the growth rate of the field in this regime 
increases as Rm 1 ' 2 . Regarding the nonlinear regime, the saturation level of the energy of the total magnetic field is independent of 
Rm when open boundaries are used. In the case of perfect conductor boundaries, the saturation level first increases as a function 
of Rm, but then decreases proportional to Rm -1 for Rm <; 30, indicative of catastrophic quenching. These results suggest that the 
shear-induced magnetic helicity flux is efficient in alleviating catastrophic quenching when open boundaries are used. The horizontally 
averaged mean field is still weakly decreasing as a function of Rm even for open boundaries. 
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1. Introduction 

The classical view of the effect of shear in large-scale dynamos 
is that it promotes the generation of magnetic fields, and that 
it is instrumental in exciting oscillatory solutions (e.g. Parker 
1955; Steenbeck & Krause 1969). While these aspects remain 
unchanged, a much more subtle effect of shear has been discov- 
ered in studies of the saturation mechanism of the large-scale 
dynamo. If the saturation of the dynamo occurs due to magnetic 
helicity conservation, a fully periodic or closed system experi- 
ences catastrophic quenching with the energy of the large-scale 
magnetic field decreasing in inverse proportion to the magnetic 
Reynolds number (Vainshtein & Cattaneo 119921 Cattaneo & 
Hughes 1996; Brandenburg 2001). A way to avoid catasnophic 
quenching is to allow magnetic helicity fluxes to escape from the 
system through open boundaries, thus allowing the large-scale 
fields to saturate at equipartition strength in a dynamical time 
scale (Blackman & Field 120001 Kleeorin et al. 120001 One of 
the most promising mechanisms, introduced by Vishniac & Cho 
(2001), operates by driving a flux of magnetic helicity along the 
isocontours of shear. Forced turbulence simulations with shear 
have confirmed that the quenching of the a-effect is up to 30 
times weaker when closed boundaries are replaced by open ones 
(Brandenburg & Sandin l2004l . 

Recent numerical experiments with forced turbulence 
(B randen burg l2005l Yousef et al.|2008a 2008b; Brandenburg et 
al.|2008 ) and convection (Kapyla et al. 2008] hereafter Paper I; 
Hughes & Proctor 2009; see also Kapyla et al. 2010) with im- 
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posed large-scale shear flows have established the existence of 
a robust large-scale dynamo. The origin of this dynamo is still 
under debate (shear-current versus incoherent a-shear effects). 
This is discussed in detail elsewhere (e.g. Brandenburg et al. 
2008; Kapyla et al. 2009b). In forced turbulence simulations 
with fully periodic boundaries, that do not allow magnetic he- 
licity fluxes, large-scale magnetic fields show slow saturation on 
a diffusive timescale (Brandenburg 2001). Earlier results from 
convection also suggest that no appreciable large-scale mag- 
netic fields are seen if magnetic helicity fluxes are suppressed 
(Paper I). 

In most of the simulations of Paper I the isocontours of shear 
were vertical and thus intersecting the boundaries. However, the 
boundary conditions (be) imposed on the magnetic field on the 
vertical borders would either allow (open vertical-field be) or 
inhibit (closed perfect-conductor be) a net flux out of the do- 
main. It was shown that when the flux is suppressed, only weak 
large-scale magnetic fields occur in the saturated state at large 
Rm. Conversely, allowing a flux produced a significant large- 
scale field. Furthermore, using a similar setup as in Tobias et al. 
(2008), where the isocontours of shear are horizontal, a similar 
effect was found. Thus it would appear that the shear-induced 
flux plays a critical role in exciting a large-scale dynamo in tur- 
bulent convection. This has given rise to the impression that the 
dynamos seen in Paper I are solely due to the magnetic helic- 
ity flux and thus inherently nonlinear (Hughes & Proctor 2009 ). 
In the present paper we argue that this interpretation is incor- 
rect and that the lack of appreciable large-scale magnetic fields 
in some of our earlier cases with vanishing helicity flux is due 
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to the fact that the critical dynamo number for the excitation of 
the dynamo is larger in that case. We also note that, more re- 
cently, large-scale dynamos have been found in turbulent con- 
vection without shear (Kapyla et al. |2009a|, irrespective of the 
presence of magnetic helicity fluxes. 

The origin of the dynamos reported in Paper I was inter- 
preted in the context of turbulent dynamo theory in Kapyla et al. 
(2009b). An important feature of turbulent dynamos is that ac- 
cording to standard mean-field theory neither the turbulent trans- 
port coefficients nor the growth rate of the field should depend 
on the microscopic resistivity. The former prediction was con- 
firmed by Kapyla et al. (2009b), but the latter has not yet been 
studied in detail in direct simulations. This is one of the goals of 
the present paper. 

Another aim is to study the saturation behaviour of the large- 
scale dynamo with both open and closed boundary conditions. 
Firstly, in homogeneous systems, i.e. where the turbulent trans- 
port coefficients do not vary in space, full saturation is expected 
to occur on a slow diffusive time scale with closed bound- 
aries whereas with open boundaries the saturation is expected 
to happen on a dynamical time scale if a helicity flux is present 
(Brandenburg 2001). However, in inhomogeneous systems the 
situation is likely to be more complex (e.g. Mitra et al. 2010). 
Furthermore, simulations of helically forced isotropic turbulence 
have shown that, in the absence of shear, the saturation level of 
the magnetic field scales as B / B% q oc Pan -1 for open bound- 
aries (Fig. 2 of Brandenburg & Subramanian 2005a). When 
shear is present, the shear-induced magnetic helicity flux should 
alleviate the Rm-dependence of the saturation level. Preliminary 
results in Paper I indicate that this might indeed be true, but the 
results are not fully conclusive because other parameters, such 
as the Rayleigh number and the fluid Reynolds number changed 
when Rm was varied. In the present paper we perform a more 
detailed study where only Rm is varied. 



2. Model and methods 

2. 1. Setup and boundary conditions 

We use the same setup as in Paper I in which a small rectan- 
gular portion of a star is modelled by a box situated at colati- 
tude 9. The coordinate system is such that (x, y, z) correspond 
to (9, <fi, r) in a spherical coordinate system. The dimensions of 
the domain are (L x , L y , L z ) — (4. 4, 2)d, where d is the depth 
of the convectively unstable layer, which is also used as our unit 
length. The box is divided into three layers, an upper cooling 
layer, a convectively unstable layer, and a lower stable overshoot 
layer. Gravity points in the negative z direction, g = (0, 0, —g), 
rotation is in the positive z direction, O = (0, 0, 51), i.e. 9 = 
corresponding to the north pole, and shear is in the y direction, 
with U = (0, S x, 0), where S is the shear rate. 

The values (zi, Z2, Z3, 2:4) = (—0.85, 0, 1, 1.15)d give the 
vertical positions of the bottom of the box, the bottom and top of 
the convectively unstable layer, and the top of the box, respec- 
tively. We use a heat conductivity profile such that the associ- 
ated hydrostatic reference solution is piecewise polytropic with 
indices (mi, TO2, TO3) = (3, 1, 1). We apply a cooling term in 
the thin uppermost layer which makes this region nearly isother- 
mal and hence stably stratified. The bottom layer is also stably 
stratified, and the middle layer is convectively unstable. 

At the vertical boundaries we use stress-free boundary con- 
ditions for the velocity, 



For the magnetic field either vertical-field or perfect-conductor 
boundary conditions are used. Occasionally we refer to them 
also as open and closed, respectively. Thus, we have 



B X = By — 



Bx,Z By^ Z 



B, 



(vertical-field, VF, or open), (2) 
(perfect-conductor, PC, closed), (3) 



the former allowing magnetic helicity flux while the latter one 
does not. This also motivates the usage of the names 'open' and 
'closed' for the two boundary conditions. The j/-direction is peri- 
odic and shearing-periodic conditions are used in the x-direction 
(e.g. Wisdom & Tremaine ll9881 ). A constant temperature gradi- 
ent is maintained at the base which leads to a steady influx of 
heat due to the constant heat conductivity. 

The simulations were performed with the PENCIL Cod^3, 
which uses sixth-order explicit finite differences in space and 
third-order accurate time stepping method. Resolutions of up to 
512 3 mesh points were used. 

2.2. Units, nondimensional quantities, and parameters 
Dimensionless quantities are obtained by setting 



Po = cp = po 



1 



(4) 



where g is the constant acceleration due to gravity, po is the den- 
sity at Z2, cp is the heat capacity at constant pressure, and fio 
is the vacuum permeability. The units of length, time, velocity, 
density, specific entropy, and magnetic field are then 



[x] = d, [t] = 
[s] = c P , [B] = \J dgp p Q 



% , [p] = Po , 



(5) 



The system is controlled by the Prandtl, Reynolds, and Rayleigh 
numbers, 



Pr 
Re- 

Ra 



Pm 



Xo 

^rms 



gd 4 ( 1 ds 
vxo \ c P dz 



Rm = rc* 8 = Pm Re 



rjk f 



(6) 
(7) 

(8) 



where v is the kinematic viscosity, \o = K/ (p m cp) is the ther- 
mal diffusivity with p m the density at the middle of the convec- 
tively unstable layer at z m = \{zj, — Z2), r\ is the magnetic diffu- 
sivity, u lms the root mean square velocity, and fcf = 2-k jd is an 
estimate of the wavenumber of the energy-carrying eddies. The 
entropy gradient, measured at z m in the non-convecting initial 
state, is given by 



1 ds 

cp dz 



V- Vad 

Hp 



(9) 



with V ad = 1 - 1/7 and V = (<91nT/cHnp) Zm , where 
7 = cp/cy = 5/3 is the ratio of the heat capacities at constant 
pressure and volume, respectively, and Hp is the pressure scale 
height, also at z = z m . The differentials of logarithmic temper- 
ature and entropy are related via ds/cp = dlnT — V a ddlnp, 
and d \np = d In T + d In p. 

The degree of stratification is controlled by the parameter 



Co 



(7 ~ l)eo 
gd 



(10) 



U v 



u z = o. 



(1) 



http : / /pencil-code . googlecode . com/ 
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Table 1. Summary of the runs with varying Pm. 



Run 


grid 


Ma 


Pm 


Rm 


A 






-Srms 


B x 


By 


B 


be 


VF1 


128 3 


0.049 


0.05 


0.78- 


-0.027 


-0.070 


1.63 


- 


- 


- 


- 


VF 


VF2 


128 3 


0.034 


0.10 


1.1 


0.013 


-0.073 


1.58 


1.04 


0.14 


0.99 


1.00 


VF 


VF3 


128 3 


0.041 


0.17 


2.2 


0.030 


-0.067 


1.58 


2.04 


0.19 


1.98 


1.99 


VF 


VF4 


128 3 


0.046 


0.25 


3.7 


0.036 


-0.066 


1.60 


2.28 


0.22 


2.18 


2.18 


VF 


VF5 


128 3 


0.046 


0.50 


7.4 


0.042 


-0.066 


1.60 


2.02 


0.19 


1.85 


1.85 


VF 


VF6 


128 3 


0.044 


1.00 


14 


0.046 


-0.068 


1.63 


1.94 


0.22 


1.65 


1.66 


VF 


VF7 


128 3 


0.042 


2.50 


34 


0.038 


-0.068 


1.60 


2.06 


0.22 


1.62 


1.63 


VF 


VF8 


256 3 


0.040 


5.00 


63 


0.032 


-0.069 


1.63 


2.13 


0.20 


1.57 


1.58 


VF 


VF9 


256 3 


0.039 


10.00 


122 


0.047 


-0.062 


1.64 


2.09 


0.20 


1.31 


1.32 


VF 


VF10 


512 3 


0.035 


20.00 


222 


0.071 


-0.067 


1.68 


2.14 


0.21 


1.09 


1.11 


VF 


PCI 


128 3 


0.045 


0.25 


3.6 


-0.019 


-0.065 


1.61 


— 


— 


— 


— 


PC 


PC2 


128 3 


0.034 


0.50 


5.4 


0.011 


-0.072 


1.57 


1.08 


0.14 


0.92 


0.93 


PC 


PC3 


128 3 


0.035 


0.67 


7.5 


0.016 


-0.072 


1.58 


1.67 


0.20 


1.49 


1.50 


PC 


PC4 


128 3 


0.040 


1.00 


13 


0.021 


-0.076 


1.58 


2.44 


0.26 


2.23 


2.23 


PC 


PC5 


128 3 


0.040 


2.50 


32 


0.028 


-0.072 


1.58 


3.22 


0.31 


2.90 


2.90 


PC 


PC6 


256 3 


0.038 


5.00 


61 


0.027 


-0.062 


1.62 


3.05 


0.33 


2.55 


2.56 


PC 


PC7 


256 3 


0.035 


10.00 


112 


0.049 


-0.065 


1.63 


2.29 


0.22 


1.53 


1.55 


PC 


PC8 


384 3 


0.037 


14.29 


168 


0.057 


-0.056 


1.65 


1.83 


0.17 


0.95 


0.96 


PC 


PC9 


512 3 


0.038 


20.00 


239 


0.064 


-0.054 


1.64 


1.62 


0.16 


0.74 


0.76 


PC 



Notes. A, Cf , and kf 1 ^ are given for the kinematic stage whereas all the other numbers are given for the saturated state of the dynamo. The last 
column specifies the vertical boundary condition for the magnetic field. In all runs, Co « 0.8, Sh « —0.4, Re « 12, and Pr « 0.7. The adiabatic 
sound speed squared varies in the range cl/(gd) = 0.33 — 1.38 in the domain. Six of the runs are discussed in particular detail, so we have marked 
them here in boldface and have also highlighted the corresponding values of Rm. 



Table 2. Summary of the runs with varying Co and Sh. 



Run 


grid 


Co 


A 




k^/ki 


be 


VF11 


128 3 


0.07 


0.000 


-0.014 


1.29 


VF 


VF12 


128 3 


0.16 


0.006 


-0.022 


1.40 


VF 


VF13 128 3 


0.34 


0.023 


-0.041 


1.52 


VF 


VF14 


128 3 


0.63 


0.038 


-0.068 


1.60 


VF 


VF15 


128 3 


1.36 


0.054 


-0.097 


1.76 


VF 


PC9 


128 3 


0.07- 


-0.008 


-0.012 


1.30 


PC 


PC10 


128 3 


0.16 


0.003 


-0.022 


1.41 


PC 


PC11 


128 3 


0.38 


0.014 


-0.042 


1.50 


PC 


PC12 


128 3 


0.75 


0.028 


-0.072 


1.58 


PC 


PC13 


128 3 


1.55 


0.034 


-0.098 


1.74 


PC 



Notes. Here Ma w 0.045, Rm « 35, Pm = 2.5, and Sh = -±Co 
(i.e. S = — fi) in all runs. All numbers are given for the kinematic state. 



where cq is the internal energy at z = z^. The internal energy 
is related to the temperature via e = cyT. We use £ = 1/3 in 
all models, which yields a density contrast of about 27 over the 
whole domain. 

The strengths of rotation and shear are quantified by the 
Coriolis and shear numbers 

2P. S 
Co= -, Sh= (11) 

where (w rms fcf) _1 gives an estimate of the convective turnover 
time. 



2.3. Characterization of the simulation data 

Our runs can be characterized using the following quantities. We 
express the rms velocity in the unit of (gd) 1 / 2 , which is related 
to the Mach number, and thus define 

Ma = u^/igd) 1 ' 2 . (12) 
The normalized growth rate of magnetic field is expressed as 

A = A/(u rms fc f ). (13) 

The typical wavenumber of the energy-carrying eddies and the 
relative kinetic helicity are given by 

fc f (w) = w rms /?t rms , ef = (u> ■ u)/((u lms )(uj lms )), (14) 

where u> = V x u is the vorticity and the angular brackets denote 
volume averaging. The normalized rms magnetic field strength 
is given by 

^rms -^rms/^eqi (15) 

where B cq — y Jiq (pU 2 ) is the equipartition value of the field. 
Corresponding values for the horizontally averaged (denoted by 
an overbar) horizontal components of the field are given by 

% = {B 2 x ) l / 2 /B c ^ \ = {B 2 y ) 1/2 /B cq , (16) 

and B = (B x + By) 1 ' 2 , is referred to as the total mean field. 

Error bars are estimated by dividing the time series into three 
equally long parts and computing a mean value individually for 
each of these. The largest deviation from the average over the 
full time series is taken to represent the error. 
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3. Results 

In order to study the effects of different boundary conditions 
and magnetic helicity fluxes on the excitation and saturation 
of convective dynamos, we perform a number of simulations 
with perfect-conductor and vertical-field boundary conditions. 
We follow two lines in parameter space where we either vary 
Pm and Rm, keeping all other parameters fixed, or we fix both 
Reynolds numbers and vary Co and Sh. The runs are listed in 
Tables Q]and|2] 

3.1. Dynamo excitation 

We begin by investigating the linear regime where the magnetic 
field is still weak. In that case the values of it xms are not yet 
quenched by the field, and so the resulting values of Rm are 
somewhat larger than in the saturated case. Turbulent dynamo 
theory predicts that large-scale dynamo action is possible if the 
magnetic Reynolds number and a relevant dynamo number ex- 
ceed critical values that depend on the details of the system in 
question (e.g. Moff att [19751 Krause & Radler [19801 Rudiger 
& Hollerbach 2004). Furthermore, boundary conditions also in- 
fluence the excitation of large-scale dynamos (e.g. Choudhuri 
1984; Jouve et al. 2008). In order to study these effects, we de- 
termine the critical Rm and the critical dynamo number, here im- 
plicitly described by the strengths of rotation and shear, needed 
to excite the dynamo with both open and closed boundary con- 
ditions. 

3.1.1. Dependence on Rm 

From the perspective of mean-field dynamo theory a, large-scale 
dynamo should operate if Rm 1 . We start the simulations with 
a random small-scale magnetic field of the order of 10~ 4 £? cq and 
monitor the time evolution of the total magnetic field quantified 
by B rms . We indeed find that in runs with vertical-field bound- 
ary conditions the dynamo is excited for Rm i; 1, whereas in 
the case of perfectly conducting boundaries, the critical value is 
about five times higher (see Table[T]and Fig. [TJ. 

As a consistency check we use a one-dimensional mean-field 
model written in terms of the vector potential 

Ai = -U jti Aj + aijBj - (rjij + //<>,, )/>,.../ j (17) 

where Ai denotes a time derivative and Uj^Aj — (SA y ,0, 0) 
is the shear term. The mean magnetic field is given by B = 
(— A y , A x , 0), the mean current density is given by /ioJ s : = 

— A i , and primes denote z-derivatives. The coefficients and 
rjij, describing the turbulent transport coefficients relevant for 
z-dependent large-scale fields, are taken from a test field calcu- 
lation (Run D4 of Kapyla et al. 2009b; see their Fig. 13) which 
is the same as our run VF7, but with Pm = 5 instead of 2.5 
used here. We then use either VF or PC boundary conditions, 
defined by Eq. (0, and determine the marginal value of 77 that 
allows dynamo action. This is quantified by a Reynolds number 
Rm c = M rms / (rjcki), where it rms is taken from the same run as 
the turbulent transport coefficients and r/ c is the marginal value 
of the microscopic magnetic diffusivity. 

We find a marginal value Rm c « 1.0 for the vertical-field 
boundary conditions, whereas in the case of perfect-conductor 
boundaries Rm c « 4.1 is obtained. These values are very close 
to those obtained from direct simulations, although one must 
bear in mind that the Rm-dependence of ay and rjij at small 



0.10 



0.01 





Vertical Field 




■ 

/ 


/ 


Rm 1/a 






Rm// Perfect 


Conductor 








1 
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Rm 
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Fig. 1. Growth rate A of the magnetic field as a function of Rm 
from the runs listed in Table Q] Solid (dashed) line with closed 
(open) symbols shows the results for the VF (PC) boundary con- 
ditions. Dotted lines proportional to Rm and Rm ' are shown 
for reference. 

values of Rm (Kapyla et al. 2009b) is not taken into account 
in our mean-field modelling. However, the qualitative result that 
the dynamo is harder to excite when perfect-conductor bound- 
aries are used is consistent with the direct simulations. 

We measure the growth rate, A, of the rms magnetic field as 

A= <dlnfl rma /di) t , (18) 

where (. . .) t denotes a time average over the kinematic range 
where B rms grows exponentially. We find that, for Reynolds 
numbers close to the marginal value, A increases roughly pro- 
portional to Rm (see Fig.[T} in accordance with theoretical con- 
siderations. However, soon this increase ceases and A levels off 
to a roughly constant value for intermediate values of Rm. For 
Rm <; 60 the growth rate is consistent with A oc Rm 1 ' 2 . This 
behaviour can be understood as follows: for small Rm the large- 
scale dynamo is excited, but the growth rate is still affected by 
the relatively large value of the microscopic magnetic diffusiv- 
ity 77. In the present system the small-scale dynamo is excited 
when Rm i; 30, but this contribution begins to dominate over 
the growth rate of the large-scale field only for Rm J; 60 where 
our data for A is consistent with an Rm ' scaling. This is in ac- 
cordance with findings of Schekochihin et al. (2004) and Haugen 
et al. (2004) for forced isotropic turbulence. 

3.1.2. Dependence on shear and rotation 

In addition to the dependence on the magnetic Reynolds number, 
we study the influence of shear on the excitation of the dynamo. 
We use a setup where Rm « 35 in the case Co = Sh = 0, 
in which case the small-scale dynamo is marginal. We then in- 
troduce shear and rotation into the system and keep the quantity 
q = — S/ft = 1 fixed in all cases. Our results for the growth 
rate of the magnetic field in the range Co ss 0.07. . . 1.5 are 
shown in Fig. [2] We find that for a given Co the growth rate 
is always greater in the simulations with open boundaries (see 
also Table|2]). The case Co ss 0.07 (Run VF11) is very close to 
marginal in the case of open boundaries, but clearly subcritical 
when perfectly conducting boundaries are used (Run PC9). The 
critical value of Co in the perfect-conductor case seems to be 
roughly twice as large as the corresponding value in the vertical- 
field runs. To be more precise, we estimate the relevant dynamo 
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£ 0.02 
< 
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0.1 



1.0 
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Fig. 2. Growth rate A of the magnetic field as a function of Co 
with open (VF, solid line) and closed (PC, dashed line) boundary 
conditions for the runs listed in Table [2] 
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Fig. 3. Dynamo numbers according to Eq. d23l as functions of 
rotation for the runs listed in Table [2] The solid (dashed) curve 
shows the results for vertical-field (perfect-conductor) boundary 
conditions. The horizontal dash-dotted and dash-triple-dotted 
lines indicate the critical dynamo number in the case of open 
and closed boundaries, respectively. Power law proportional to 
Co 2 is shown for reference. 



numbers as 

m 

C a = 

Cs = 



iltki 
S 



kx ' 



3Sh 



(19) 
(20) 



where fci — n/L z is the wavenumber of the leading decay mode 
for both boundary conditions. Here we approximate the a effect 
and turbulent diffusivity by 



■\t c {u)-u), 



±T U 2 

3 ' c "-rms ■ 



Furthermore, we assume the Strouhal number, 

St = T c lt rms fcf , 



(21) 



(22) 



of unity (as found by Brandenburg & Subramanian 2005b) so 
that r c = (unnsfcf) -1 . The fractional helicity is defined via ef = 



(uj-u)/(u lms uj rms ), where uj t 
the dynamo number as 

D = C a C s . 



- k {uj) n 
ms "'f u y 



. We can now define 



(23) 



For the marginal open-field Run VF1 1 we find a critical dynamo 



number D 



(VF) 
crit 



0.1. The dynamo numbers for the perfect con- 
ductor Runs PC9 and PC 10 are 0.7 and 3.5, respectively. Linear 
interpolation of the magnetic field growth rates give a critical 



value of D 



(PC) 
crit 



0.3 which is roughly three times greater than 



in the vertical-field case (see Fig. [3). We find that the dynamo 
number is approximately proportional to Co 2 in the range of pa- 
rameters studied here. This dependence can be understood qual- 
itatively as follows: the a effect, and hence C a , increases ap- 
proximately proportional to Co in the range of parameters stud- 
ied here (see Kapyla et al. 2009b). Furthermore, as we keep the 
ratio S/ £1 constant, shear and Cs also increase proportionally to 
Co. 

In the preceding simplified analysis we have estimated a 
and rjt using volume averages that are not ideally suited for 
the present case (see, e.g. Kapyla et al. |2009bl >. Furthermore, 
we have assumed here that the dynamo is of pure a-shear type 
and have neglected all other possible induction effects (shear- 
current, fix J, and incoherent a-shear effects) that can assist in 
dynamo action (Kapyla et al. 2009b). This means that the quan- 
titative results should be taken with some amount of caution. 

In Paper I, the only case where open and closed boundaries 
were compared was the one where only shear was present, i.e. 
Sh 0, whereas Co = (Run C with vertical-field bound- 
ary conditions versus Run C with perfect-conductor boundary 
conditions). In Paper I we also showed that, in the absence of 
rotation, the growth rate of the magnetic field is reduced in com- 
parison to the case where both shear and rotation are present. 
Furthermore, our present results suggest that the large-scale dy- 
namo is harder to excite in the case of perfect-conductor bound- 
aries. Thus, a reasonable explanation to the slow growth of the 
large-scale field in Run C of Paper I is that the system is close 
to being marginally excited, and not, as suggested by Hughes 
& Proctor (2009), the result of some nonlinear effect arising 
from magnetic helicity flux out of the system. A straightforward 
quantitative analysis in terms of a classical a-shear dynamo de- 
scribed above is, however, not useful in the non-rotating case. 
This is shown by the results of Kapyla et al. (2009b): the main 
contribution to dynamo action in this case is due to the incoher- 
ent a-shear process which relies on the fluctuations of a (e.g. 
Vishnia c & B randenburg [19971 Sokolov [19971 Silant'ev[2000; 
Proctor l2007t . 

3.2. Nonlinear regime 

We find that in all cases where a dynamo is excited, large-scale 
magnetic fields are also generated. In the low-Rm runs the large- 
scale contribution is substantial, i.e. B/B ims is close to unity al- 
ready in the kinematic regime. In the cases where also the fluctu- 
ation dynamo is excited, a large-scale pattern is discernible only 
once the dynamo has reached saturation. This feature is com- 
mon to all known large-scale dynamos and was already seen in 
simulations of Brandenburg (2001). Figure [4] shows visualiza- 
tions of the streamwise component of the field in the saturated 
state for three Reynolds numbers, representing low, intermedi- 
ate, and high Rm, with open and closed boundaries (see Table 
Q~|i. The runs depicted in Fig. [4] are indicated in boldface in the 
table. Although an increasing amount of small-scale features is 
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4,0 



Fig. 4. Magnetic field component B y from runs with vertical-field (left column) and perfect-conductor (right column) boundary 
conditions. Three different cases with low (top), intermediate (middle) and high (bottom) Rm are shown. The sides of the boxes 
show the By field at the periphery of the domain whereas the top and bottom panels show the field at z — z$ and z — z%, respectively. 



seen with increasing Rm, a large-scale pattern is clearly visible 
in all cases. 

Unlike the case with just rotation and no shear, where the 
mean field shows variations in the horizontal directions (Kapyla 
et al. 2009a), here the k x — k y — mode dominates the large- 
scale field and thus horizontal averages are suitable. Space-time 
diagrams of the horizontally averaged magnetic field compo- 
nents, B x and B y , are shown in Fig. [5] for the same runs as in 
Fig. |4] We find that the large-scale field is non-oscillatory in all 
cases, which is in agreement with earlier results (Paper I; Hughes 
& Proctor 2009; Kapyla et al. 120101 1. In most cases, regardless 
of the boundary conditions, the field changes sign near the base 
of the convectively unstable layer, with the exception of the runs 
with vertical-field boundary conditions (VF runs) with the low- 
est values of Rm. (The solutions are invariant under sign rever- 
sal, B — > B, and both realizations of the large-scale field are 
found in Figs. [4] and [5] depending just on the initial conditions.) 



In the perfect-conductor runs (PC runs), the field near the top of 
the domain has a different sign than in the bulk of the convection 
zone, whereas in the VF runs such behaviour is not seen. We find 
that in the PC runs the layer of oppositely directed field becomes 
progressively thinner as Rm is increased, leading to strong gra- 
dients near the boundary (see Fig . O . This can possibly explain 
the numerical problems encountered in the Rm « 239 simu- 
lation (Run PC9) with perfect-conductor boundary conditions. 
The saturated state of Run PC9 is significantly shorter than that 
of the other runs due to numerical instability that prevented run- 
ning the simulation further. 

It turns out that the emergence of the large-scale field oc- 
curs progressively earlier as Rm is increased. This can be seen 
most clearly in Fig. [5] which shows that the large-scale field has 
reached saturation in less than 200 turnover times for Rm = 222 
(Run VF10), while for Runs VF7 and VF4 with Rm = 34 and 
3.7, the saturation times exceed 300 and 400 turnover times, re- 
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Fig. 5. Horizontally averaged magnetic field components B x and B y from the same six runs as in Fig. |4] Those runs and the 
corresponding values of Rm are also highlighted in Table Q] in bold face. The dotted white lines at z = and z = d denote the 
bottom and top of the convectively unstable layer, respectively. The runs on the left are for open boundaries and those on the right 
are for closed ones. The magnetic Reynolds number increases downwards. 



spectively. A similar trend is seen also in the PC runs. In an ear- 
lier study we found that the mean values of the turbulent trans- 
port coefficients relevant for the generation of large-scale mag- 
netic fields, ctyy and rj t , remain constant within the errors given 
that Rm J; 8; see Fig. 5 of Kapyla et al. (2009b). However, as 
the magnetic Reynolds number is increased, the fluctuations of 
a tend to increase as well. Such fluctuations can contribute to 
the incoherent a-shear process which can possibly explain the 
faster saturation of the large-scale field when Rm is increased. 



The saturation level of the total and mean magnetic fields 
for runs with open and closed boundary conditions are shown 
in Fig. [6] We find that the total magnetic field in the runs with 
open boundary conditions is roughly consistent with an Rm- 
independent value. In the perfect conductor case the total mag- 
netic energy first increases up to Rm w 30 and then decreases 
proportionally to Rm -1 for Rm ^ 60. Mean-field models tak- 
ing into account magnetic helicity evolution can also produce a 
maximum for the saturation level at some intermediate Rm (e.g. 
Brandenburg et al. 2007). We note that an Rm -1 dependence 
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for the mean field energy is indicative of catastrophic quench- 
ing. This would indeed be the expected result for closed bound- 
aries. However, in our case only the total field energy shows the 
Rm" 1 behaviour whereas the mean field exhibits a much steeper 
(at least Rm ) dependence. The explanation for such a steep 
trend is as yet unclear. The data for the mean magnetic field in 
the case of open boundaries also show a weak decreasing trend 
consistent with a power law Rm as opposed to the expec- 
tation that the saturation level is independent of Rm. However, 
the unexpected behaviour of the saturation level could simply 
be related to the fact that in mean-field models true asymptotic 
behaviour may only commence at much larger values of Rm 
(Brandenburg et al. l2009b . 

Our simulations with the highest magnetic Reynolds num- 
bers and closed boundaries apparently do not show a slow satu- 
ration behaviour; e.g. as in Fig.|7]where the mean magnetic field 
and a saturation predictor (Brandenburg !2Q0Tb proportional to 
1 — e~ 2r,k ( t_t = at ) are shown for Runs PC6-PC8. Here we use 
k = 2ir/L z and t sat is the time at which the small-scale dy- 
namo has saturated. In the runs with intermediate Rm (PC6 and 
PC7) the saturation predictor is in fairly good agreement with 
the simulation results, whereas for Run PC8 this is no longer 
the case. This might be caused by additional contributions to B 
whose effective value of k is larger. In fact, earlier simulations 
of forced turbulence with perfect conductor boundary conditions 
(Brandenburg & Dobler 2002 ) showed that the final configura- 
tion of the mean field can be established in steps, but the time 
between different steps can still be resistively long. This could 
explain why, in the highest-Rm simulations, the saturation level 
of the mean magnetic field is lower than expected. Whether this 
is also the case in the present simulations remains an open ques- 
tion. 



4. Conclusions 

We have studied the effects of magnetic boundary conditions on 
the excitation and saturation of large-scale dynamos driven by 
turbulent convection, shear, and rotation by means of numerical 
simulations. We find that the critical magnetic Reynolds number 
is greater (Rm » 5) in the case of closed boundaries (perfect- 
conductor) compared to simulations with open (vertical-field) 
boundaries where the large-scale dynamo is excited already for 
Rm w 1. A similar result is obtained from mean-field mod- 
els operating in the same parameter regime. Furthermore, the 
relevant dynamo number is roughly three times greater with 
closed boundaries. These effects are likely to explain the weak 
large-scale dynamo action seen in the perfect-conductor runs in 
Paper I. 

The measured growth rate of the magnetic field is indepen- 
dent of the microscopic resistivity when Rm is sufficiently above 
the critical value. This is manifested by the approximately con- 
stant growth rate in the intermediate Rm range in Fig. Q] For 
Rm 30, the small-scale dynamo is excited, and for Rm 60 
it becomes dominant. The growth rate of the magnetic field is 
then consistent with A oc Rm 1 / 2 scaling, which is in accordance 
with the results of Schekochihin et al. (2004) and Haugen et al. 
(2004). The fact that the qualitative behaviour of the growth rate 
of the magnetic field is similar for both boundary conditions sug- 
gests that the origin of the large-scale dynamo is not likely to be 
a process that is essentially nonlinear, as suggested by Hughes 
& Proctor (2009), but that it can be understood within the frame- 
work of classical kinematic mean-field theory. 
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Fig. 6. Upper panel: rms-value of the total magnetic field in the 
saturated regime for runs with perfect-conductor (solid line) and 
vertical-field (dashed) boundaries. Lower panel: same as above 
but for the horizontally averaged mean magnetic field B. In both 
panels power laws fitting to the data at high Rm are shown for 
reference. 



In the saturated state the energy of the total magnetic field 
remains independent of Rm for open boundaries and decreases 
as Rm -1 for closed boundaries. The latter result is consistent 
with catastrophic quenching while the former result suggests that 
magnetic helicity fluxes are efficiently driven out of the system. 
On the other hand, the energy of the mean field, taken here to 
be represented by horizontal averages, decreases approximately 
as Rm" 25 for open and as Rm for closed boundaries. It is 
not yet clear why the mean fields tend to show a steeper decline 
than the total field. It is possible that this declining trend levels 
off at a higher Rm and that the magnetic Reynolds numbers in 
our simulations are still not large enough (cf. Brandenburg et al. 
2009). A similarly weak Rm-dependence has been observed in 
the cycle period of a-shear dynamos with isotropically forced 
turbulence (Kapyla & Brandenburg 2009]). 

We find that, for intermediate values of Rm, the large- 
scale magnetic field saturates on a resistive time scale when 
closed boundaries are used. However, with the largest Rm no 
clear signs of slow saturation are observed. Earlier results us- 
ing perfect-conductor boundaries have shown that the mean 
field can evolve in steps (Brandenburg & Dobler 2002; see also 
Brandenburg et al. 2007 ) which are associated with a change of 
the large-scale magnetic field configuration. We have not seen 
such a behaviour in our current simulations but the existence of 



8 



Kapyla' et al.: Open and closed boundaries in convective dynamos 




200 400 600 800 



Fig. 7. Mean magnetic field as a function of time from Runs PC6 
(solid line), and PC7 (dashed). The inset shows the same for 
Run PC8. The dotted lines show a saturation predictor according 
to Brandenburg (2001) with the microscopic values of r\ and k = 

2w/L z . 
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such events at a later stage cannot be ruled out. Another possi- 
ble explanation is that there are magnetic helicity fluxes occur- 
ring inside the domain which arise from the spatial gradients of 
magnetic helicity (e.g. Covas et al. |1998t Kleeorin et al. 2000; 
Mitra et al . 120 1 01 > . However, a quantitative study of these effects 
requires more detailed knowledge of the helicity fluxes and pos- 
sibly an anisotropic formulation of the magnetic a effect. These 
issues merit further investigation and are beyond the scope of the 
present paper. 
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